Quantum Chaos Border for Quantum Computing 
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We study a generic model of quantum computer, composed of many qubits coupled by short-range 
interaction. Above a critical interqubit coupling strength, quantum chaos sets in, leading to quantum 
ergodicity of the computer eigenstates. In this regime the noninteracting qubit structure disappears, 
the eigenstates become complex and the operability of the computer is destroyed. Despite the fact 
that the spacing between multi-qubit states drops exponentially with the number of qubits n, we 
show that the quantum chaos border decreases only linearly with n. This opens a broad parameter 
region where the efficient operation of a quantum computer remains possible. 

PACS numbers: 03.67.Lx, 05.45.Mt, 24.10.Cn 



Since the pioneering work of Feynman Jl| and mod- 
ern developments of efficient algorithms |2J and error- 
correcting codes the realization of quantum com- 
puters became a challenge of modern physics ||. Dif- 
ferent experimental realizations have been proposed, in- 
cluding ion traps || , nuclear magnetic resonance systems 
0, nuclear spins with interaction controlled electroni- 
cally > quantum dots Q , Cooper pair boxes and 
optical lattices |l2) . A key common feature of these ex- 
perimental settings is the presence of interacting qubits 
(two- level systems). Here we analyze the effect of qubit 
interaction on operability of the quantum computer. The 
interaction is required since a quantum computer needs 
to perform two-qubit logical operation such as XOR || . 
We note that such a two-qubit gate has been experimen- 
tally realized (b| . 

In an isolated system of n uncoupled qubits, the di- 
mension of the total Hilbert space Njj increases expo- 
nentially with n (Njj = 2"), while all eigenvalues of the 
Hamiltonian are included in an energy interval of size 
AE ~ nAo, where Ao is the average energy distance 
between the two states of one qubit. As a result, the 
average spacing A„ between adjacent energy levels of 
the Hamiltonian decreases exponentially with the num- 
ber of qubits (A„ ~ nAo/Njj -C A ). When a cou- 
pling J between the qubits is added (J < A ), one still 
has AE ~ nAo, Njj is unchanged, and the above esti- 
mate for A„ still holds. This general result for A„ is 
related to the exponentially large size Njj of the Hilbert 
space, which is one of the main reasons of the striking 
efficiency of quantum computing 0,^]. It implies that 
dense highly excited states are needed for the computa- 
tion. However, when performing the computation one 
wants to operate with noninteracting multi-qubit states 
\ipi >= \ai,...,On > where = 0,1 marks the polar- 
ization of each individual qubit. These quantum regis- 
ter states should remain well-defined in the presence of 
interqubit coupling even if multi-qubit levels are expo- 
nentially dense. Therefore the mixing of noninteracting 
multi-qubit states induced by the interaction is crucial for 



the computer operability. In the field of quantum chaos 
|l4] , ^5| it is known that noninteracting states will be even- 
tually mixed by the interaction and quantum ergodicity 
will set in: each quantum computer eigenstate will be 
composed of a large number of noninteracting multi-qubit 
states \ipi > and the original quantum register states will 
be washed out. At first glance one would expect that 
such mixing happens when the coupling between qubits 
becomes comparable to the multi-qubit spacing A„. In 
such a case, the creation of quantum computers compet- 
itive with classical ones would be rather difficult: since 
hundreds of qubits are necessary, this would lead to ab- 
surdly strict restrictions on coupling strength. Indeed, 
for n — 1000, the minimum number of qubits for which 
Shor's algorithm becomes useful || , the multi-qubit spac- 
ing becomes A„ ~ 10 3 x 2~ lo3 A ~ 10~ 298 K, where we 
used Ao ~ 1 K that corresponds to the typical one-qubit 
spacing in the experimental proposals K|. It is clear 
that the residual interaction J between qubits in any ex- 
perimental realization of the quantum computer will be 
larger than this. For example, in the proposal 0, the 
increase of effective electron mass by a factor of two, in- 
duced by the electrostatic gate potential, means that the 
spin-spin interaction is changed from J ~ Ao ~ 1 K (cor- 
responding to a distance between donors of 200 A and an 
effective Bohr radius of 30 A in Eq.2 of Q) to the residual 
interaction J ~ 10~ 5 K » A„. 

However the problem is not so simple, since the in- 
teraction is always of two-body nature and not all of 
the multi-qubit states are directly coupled. Actually the 
number of states directly coupled to such a quantum reg- 
ister state \ipi > increases not faster than quadratically 
with n. A similar problem appears in other physical 
many-body interacting systems such as nuclei, complex 
atoms, quantum dots and quantum spin glasses [[[(i] po| . 
It was realized that sufficiently strong interaction leads to 
quantum chaos and internal (dynamical) thcrmalization, 
where the eigenstates properties follow the predictions of 
Random Matrix Theory (RMT) The quantum 

chaos border for this dynamical thermalization has been 
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established only recently and it has been shown that the 
relevant coupling strength should be larger than the en- 
ergy spacing between directly coupled states A c Jl7| , pot . 
Since A c drops algebraically with n, it is exponentially 
larger than A„ ~ ?i2~ n Ao, and therefore a relatively 
large coupling strength is required for the emergence of 
quantum chaos and ergodicity. A similar border for in- 
teracting qubit systems would allow a reasonable regime 
of operability for quantum computers. 

To investigate the emergence of quantum chaos in 
quantum computers, we chose a model of n qubits 
on a two-dimensional lattice with nearest-neighbour in- 
terqubit coupling. The Hamiltonian reads: 

i i<j 

where the at are the Pauli matrices for the qubit i and 
the second sum runs over nearest-neighbour qubit pairs 
with periodic boundary conditions applied. The energy 
spacing between the two states of a qubit is represented 
by randomly and uniformly distributed in the interval 
[A — 5/2, Aq + 5/2]. The parameter 5 gives the width 
of the distribution near the average value Ao and varies 
from to Ao. Here 1^ can be viewed as the splitting 
of nuclear spin levels in a local magnetic field, as it is 
discussed in the experimental proposals |^]|]] . The differ- 
ent values of I\ are needed to prepare a specific initial 
state by electromagnetic pulses in nuclear magnetic res- 
onance. In this case the couplings Jy will represent the 
hyperfine interaction between the spins, which is needed 
to build the quantum computer. Different physical mech- 
anisms can generate these couplings, such as spin exci- 
ton exchange j8|J|, dipole-dipole interaction, etc... For 
generality we chose randomly distributed in the in- 
terval [— J, J]. The Hamiltonian (|]) can be considered as 
a generic quantum computer model, which catches the 
main physics of different experimental proposals. For 
example a similar Hamiltonian appears in a quantum 
computer based on optical lattices |TJ,|2l]]. We restrict 
ourselves to the case of static couplings which are always 
present as a residual interaction and are much larger than 
the multi-qubit spacing A„ even for moderate values of 
n. In a sense (|l|) describes the hardware of the computer, 
while gates operation in time requires additional studies, 
which are possible only if the properties of the hardware 
are well understood. 

As is well known in the field of quantum chaos, 
the transition to ergodic eigenstates is reflected in 
the level spacing statistics P(s), which goes from the 
Poisson distribution -Pp(s) = exp(— s) for nonergodic 
states to the Wigner-Dyson (WD) distribution Pw{s) = 
(ns/2) exp(— 7rs 2 /4), corresponding to RMT, for ergodic 
states. Here s is the nearest level spacing measured in 
units of average spacing and P(s) is the probability to 
find two adjacent levels whose spacing is in [s, s + ds]. 
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FIG. 1. Transition from Poisson to WD statistics in 
the model (Q for the states in the middle of the 
energy band (±6.25% around the center) for n—12 : 
J/Aq — 0.02, r\ — 1.003 (dashed line histogram); 
J/A = 0.48, r] = 0.049 (full line histogram). Full curves 
show P P {s) and Pw(s); N s > 2.5 x 10 4 , N D = 100, 8 = A . 




FIG. 2. Dependence of r/ on the rescaled coupling strength 
J/J c for the states in the middle of the energy band for 
n = 6(*),9(o), 12(triangles), 15(squares); 5 = Ao. The upper 
insert shows log(J c /Ao) (diamonds) and log(J cs /Ao) (trian- 
gles) versus log(n); the variation of the scaled multi-qubit 
spacing A„/Ao with log(n) is shown for comparison (+). 
Dashed line gives the theoretical formula J c = CAo/n with 
C — 3.16; the solid line is J cs = 0.41Ao/n. The lower insert 
shows log(J cs /Ao) versus log(5/Ao) for n = 6(*),9 (o), 12 
(triangles); straight lines have slope 1. 

The majority of our data are displayed for the middle 
of the energy spectrum, where the transition starts, and 
which therefore sets the limit of operability of the quan- 
tum computer. The model (|l]) has two symmetry classes 
characterized by an odd or even number of qubits up, and 
the data are given for one symmetry class. In order to 
reduce statistical fluctuations, we use 5 < < 4 x 10 4 
random realizations of Tj and Jy , as is done usually in 
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RMT [[L5| . Eigenvalues and eigenvectors are computed by 
exact diagonalization of the Hamiltonian matrix (Q) for 
each realization. In this way the total number of spacings 
is 10 4 < N s < 1.6 x 10 5 (N s oc N D N H ). An example of 
the transition in the spectral statistics is shown in Fig.l. 

To analyze the evolution of P(s) with the coupling J, 
it is convenient to use the parameter r\ = Jq°{P{s) — 
P w {s))ds/ Jo°{Pp{s) - P w {s))ds, where s = 0.4729... 
is the intersection point of Pp(s) and P\y(s). In this way 
Pp(s) corresponds to 77 = 1, and P\v{ s ) to r? = 0. As 
is usual in the field of quantum chaos, the variation of 77 
characterizes the evolution of P(s) [^0|. The variation of 
77 with respect to J/A is presented in Fig. 2 for 5 = A 
showing that indeed rj drops from 1 to with increas- 
ing coupling strength. The transition appears to become 
sharper for larger system sizes. The typical J c value near 
which the transition takes place corresponds to interme- 
diate values of 77. We chose the condition 77(J C ) = 0.3. 
The dependence of J c on n is given in the Fig. 2. In anal- 
ogy with other many-body systems discussed in |p7| , po| , 
we expect that J c » A c ss CAo/n, where C is some 
numerical constant. Indeed, one multi-qubit state is cou- 
pled to 2n other states in an energy interval of order 6A0. 
This theoretical estimate is in agreement with the data 
of Fig. 2, with C ss 3. We stress that this critical cou- 
pling is exponentially larger than the multi-qubit level 
spacing A„ ~ n2~"Ao, as is shown on Fig. 2. For the 
case 5 <C Ao, the total spectrum at J = is composed of 
n bands with interband distance 2Ao and a bandwidth 
of y/nS. Within one band, one multi-qubit state is cou- 
pled to about n states in an energy interval of 26, so that 
J c ~ A c ~ 6 /n. This quantum chaos border is still much 
bigger than A„ ~ Vnd/(N H /n) ~ n 3 / 2 2- n 6. 

The transition in the level statistics reflects the drastic 
change in the multi-qubit structure of the eigenstates of 
(|I|). Indeed, Fig. 3 shows that for J < J c one eigenstate 
is formed only by one or few noninteracting states \tj)i >, 
while for J > J c a huge number of them are required. 
In the latter case, the computer eigenstates become a 
random mixture of quantum register states \tpi >, making 
rather difficult to perform computation. 

To study this drastic change in the structure of eigen- 
states, it is convenient to use the quantum eigenstate 
entropy S q , defined by: S q = — ^ ji Wi\og 2 Wi, where 
Wi is the quantum probability to find the noninteract- 
ing multi-qubit state \ipi > in the eigenstate \<f> > of (Q) 
(Wi = I < ij)i\<f) > | 2 ). In this way S q = if \<j> > is one 
noninteracting state (J = Q), S q = 1 if \<f> > is equally 
composed of two \tpi >, and the maximal value is S q = n 
if all 2" states contribute equally to \<f> >. The variation 
of the average quantum entropy with J is shown on Fig. 4 
for S = Ao. It shows that S q grows with J and the tran- 
sition to ergodic states with large S q takes place in the 
vicinity of J c . In addition these data show that the crit- 
ical coupling J cs at which S q — 1 is J cs w 0.13J C . The 



ratio J cs / Jc stays within 15% of the average value when 
n changes from 6 to 15, while the ratio A„/ J c varies from 
Ito3xl0 -3 (see upper insert of Fig. 2). The dependence 
of J cs on 6 is shown on the lower insert of Fig. 2; it clearly 
shows the linear decrease of J cs with S and can be well 
described by J cs = OAS/n. Naturally, the quantum chaos 
border drops to zero with 6 due to the quasidegeneracy 
inside the energy bands at J = 0. 




FIG. 3. Two quantum computer eigenstates of model 
(0) in the basis of noninteracting multi-qubit states, i.e. 
Wi — I < i/h\4> > 1 2 as a function of noninteracting multi-qubit 
energy Ei for n — 12 and 8 — Ao with J c /Ao = 0.273 (see 
text): (a) J/A = 0.02; (b) J/A = 0.48. 




FIG. 4. Dependence of the quantum eigenstate entropy S q 
on J/J c for S — Ao and n = 6(*), 9 (o), 12 (triangles), 15 
(squares); 10 4 < Ns < 1.6 xlO 5 . Insert shows the dependence 
on larger scale. 

We note that for n = 1000 and S = Ao = 1 K, only two 
multi-qubit states will be mixed at J cs ~ 0.4Ao/ti w 0.4 
mK. This critical coupling is much larger than the multi- 
qubit level spacing A„ ~ 10 -298 K. Even if the quan- 
tum border J cs corresponds to a relatively low coupling 
strength it seems reasonable that the residual interac- 
tion between qubits can stay below this threshold with 
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current technologies (but not below A n ). 




FIG. 5. The quantum computer melting induced by the 
coupling between qubits. Color represents the level of quan- 
tum eigenstate entropy S q , from bright red (S q « 11) to blue 
(S q — 0). Horizontal axis is the energy of the computer eigen- 
states counted from the ground state to the maximal energy 
(« 2nAo). Vertical axis is the value of J/Ao, varying from 
to 0.5. Here n — 12, 8 — Ao, J c /Ao = 0.273, and one random 
realization of (|l]) is chosen. 

The pictorial image of the quantum computer melting 
under the influence of the interqubit coupling J is shown 
on Fig. 5. The melting starts in the middle of the spec- 
trum (high energy) and progressively invades low-energy 
states and the whole computer, destroying its operability. 
We stress that this destruction takes place in an isolated 
system without any external decoherence process. Nev- 
ertheless the thermalization in this closed system, which 
appears because of the interqubit coupling, can mimic the 
effect of a coupling with the external world and external 
decoherence. Above the quantum chaos border an ini- 
tial register state > will spread quickly with time ]2^ ] 
over an exponential number of eigenstates of the system 
with residual interaction, destroying gates operability. 

Our studies show that the mixing of multi-qubit states 
and onset of quantum chaos induced by interqubit cou- 
pling leads to the melting of a realistic quantum com- 
puter and destruction of its operability; however, the 
quantum chaos border found for this process corre- 
sponds to a relatively strong interaction, being expo- 
nentially larger than the energy level spacing between 
multi-qubit states. We expect that below this border, 
error-correcting codes ||[| will allow to perform efficient 
quantum computing with a large number of qubits. How- 
ever we note that quantum chaos sets in very easily if the 
fluctuation amplitude S of individual qubit spacing drops 
to zero (J c oc S). 
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